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Abstract 

We report multicanonical Monte Carlo simulations of the tails of the order-parameter distribution of the two-dimensional 
Ising model for fixed boundary conditions. Clear numerical evidence for "fat" stretched exponential tails is found below the 
critical temperature, indicating the possible presence of fat tails at the critical temperature. 
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1. Introduction 

A quantity of central importance for finite-size scal- 
ing (FSS) analyses of critical phenomena is the order- 
parameter distribution p{m) [1]. Most properties of 
p{rn) at criticality are known from computer simu- 
lations [2]. Analytical information comes from field 
theoretic renormalization group calculations [3], con- 
formal field theory [4], and also a generalized classifi- 
cation theory of phase transitions [5,6]. While some 
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of the analytical predictions seem to have been cor- 
roborated by numerical simulations [6, 7], the predic- 
tions for the tails of the critical order-parameter distri- 
bution could not be confirmed. Recording of the very 
small probabilities in the tails requires special tech- 
niques such as multicanonical simulations [8]. Many 
recent studies [7, 9] have attempted this, but failed in 
establishing the true behaviour of the tails of p(m). 

In the present paper we report results of high pre- 
cision multicanonical Monte Carlo simulations for 
the two-dimensional (2D) Ising model on square lat- 
tices with fixed (i.e. all boundary spins fixed to +1) 
boundary conditions. One of the objectives is to study 
whether the order-parameter distribution obtained 
from the simulation can be considered to be asymp- 
totic with respect to system size. A secondary objec- 
tive is to study fixed (symmetry breaking) boundary 
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Fig. 1. Tail analyses for the 2D Ising model with fixed boundary conditions. 



conditions because the asymmetry should give rise to 
an asymmetry in the far tail behaviour. 



2. Analysis methods 

We consider the 2D Ising model on a square lattice 
with N = L 2 spins o-j = ±1 interacting according to 
the usual Hamiltonian 

H = -J^Uiaj, (1) 

where J > is the ferromagnetic coupling strength 
and the summation mns over a ^ nearest neigh- 

bour pairs on the lattice. The order parameter is the 
magnetization per spin, —1 < m = (1/N) X^Li a i ^ 
1. In the following we set J = 1, the Boltzmann con- 
stant to unity and denote the temperature by T. 



The probability density p(m) of the order param- 
eter depends parametrically on temperature and sys- 
tem size, p(m) — p(m; T, L). At criticality, traditional 
FSS predicts the scaling form [1] 

p(m; T = T c , L) - L^p(mL^), (2) 

where p(rh) with m = mL^l v is a universal scaling 
function. Since for the 2D Ising model v = 1 and 
(3 = 1/8, the scaling variable is thus m = mL 1 ^. 
For the analysis of the tails of p(m) we first determine 
constants A, B, C such that po(x) = Ap(B(m — C)) 
has mean zero, unit norm and unit variance. Then we 
split the peak into its left and right tails by defining 
the functions 

POl (x) = PO (Zpcak - X) for X < X pcak , (3) 

POr(x) =Po(x - Xpcak) for X > X pC ak , (4) 
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where x pca k is the position of the maximum. To ex- 
hibit stretched exponential tails we finally calculate 
the logarithmic derivatives 

^ )= dlog (-log loPOl) ; 
dlog 10 x 

where i — l,r, and plot them against y = log 10 a;. 
In this representation, tails of the form po(x) ~ 
B(x + c) 13 cxp[— A(x + c) a ] lead to a plateau at the 
value a for large x. A standard Gaussian distribution 
(l/\/27r) exp(— x 2 /2) thus approaches asymptoti- 
cally the value 2. 
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3. Results 

We performed multicanonical simulations [8] by 
flattening the probability distribution of the magneti- 
zation for system sizes up to L = 64 at three char- 
acteristic temperatures. The right tails in the high- 
temperature phase at T = 3.5 are shown in the up- 
per left plot of Fig. 1. With increasing system size 
the distribution approaches Gaussian behaviour as is 
expected from the central limit theorem for a finite 
correlation length £ = l/[- lntanh(l/T) - 2/T] = 
1.41235 .... In fact, assuming heuristically that there 
are effectively L 2 ff uncorrected spins with L e g — 
L/2£ (= 6, 12, 23 for L = 16, 32, 64), we obtain the 
very similar looking curves in the upper right plot of 
Fig. 1, explaining the finite-size effects. 

At the critical temperature T c = 2/arcsinh(l) = 
2.2691 . . ., theory predicts a right tail of the form 
x^- 1 )/ 2 exp(— x (5+1 ), where 5 = 15 is the equation 
of state exponent, leading to a plateau at 16 [7]. Our 
data shown in the lower left plot of Fig. 1 reveal a 
shoulder developing with increasing L. The conver- 
gence, however, is extremely slow [10]. The solid line 
is a fit using this theoretical prediction, 



Po (x) = a((x + c)/b) 7 exp(-((a: + c)/6) 16 ), 



(6) 



with a = 1.31, b = 8.59, and c = 7.79. 

For low temperatures the magnetization is close to 
unity. At T = 1.5 a true peak starts developing only 
for system sizes L > 16, and even for L = 32 only five 
data points exist to the right of the peak, cf. Fig. 2. We 
therefore found only weak evidence for Gaussian be- 
haviour with increasing L for the right tails [10]. Here 



Fig. 2. Rescaled order parameter distributions p$(x) on a logarith- 
mic scale for T = 1.5 and L = 16, 32, 64 for fixed and periodic 
boundary conditions. For periodic boundary conditions only the 
right half of the distribution is shown. 

the left tails shown in the lower right plot of Fig. 1 are 
more interesting. Near the peak (i.e. for small x) we 
see a narrow regime over which the curves approach 
Gaussian behaviour as one would expect again from 
the central limit theorem. In the intermediate range of 
x, the curves show a plateau at 0.5, corresponding to 
a fat stretched exponential ~ exp(— \/x). This rep- 
resents the well known droplet regime predicted ana- 
lytically [11]. We observed this stretched exponential 
tail in the distributions for low temperatures all the 
way up to T c . For periodic boundary conditions the 
stretched exponential tails eventually cross over into a 
flat bottom, reflecting phase coexistence on finite lat- 
tices governed by strip-like spin configurations. For 
fixed boundary conditions the same stretched expo- 
nential tail is terminated by a cutoff function. Visually 
the three different regimes can also easily be identified 
in the logarithmic plot of po{x) shown in Fig. 2. 



4. Summary 

To summarize, the order-parameter distribution at 
the critical point is found to approach its asymptotic 
universal scaling function extremely slowly. This is 
also true for periodic boundary conditions [10]. Our 
extrapolations discussed in detail in Ref. [10] indicate 
that the required system sizes (L > 10 5 ) are beyond 
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present day numerical resources even for the 2D Ising 
model. 
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